Topological magnon insulator with Dzyaloshinskii–Moriya interaction under the irradiation of light
Chen Liang
School of Physics and Electronic Engineering, Taishan University, Taian 271000, China

 

† Corresponding author. E-mail: zhuclweifang@sina.com

Project supported by the National Natural Science Foundation of China (Grant No. 61604106) and Shandong Provincial Natural Science Foundation, China (Grant No. ZR2014FL025).

Abstract

The topological magnon insulator on a honeycomb lattice with Dzyaloshinskii–Moriya interaction (DMI) is studied under the application of a circularly polarized light. At the high-frequency regime, the effective tight-binding model is obtained based on Brillouin–Wigner theory. Then, we study the corresponding Berry curvature and Chern number. In the Dirac model, the interplay between a light-induced handedness-dependent effective DMI and intrinsic DMI is discussed.

1. Introduction

A topological insulator is a material that has a bulk band gap like an ordinary insulator but whose surface contains conducting states.[1,2] The search for topological states of quantum matter is one of the hottest topics in condensed matter physics. The topologically protected edge states is not unique to electronic systems. A growing number of studies investigate the possibility of topological properties on magnons in two-dimensional (2D) honeycomb systems[312] due to their similarity to electrons in graphene.[13,14] The ferromagnetic spins on a 2D honeycomb lattice can be topologically nontrivial when a proper nearest neighbor exchange exists.[3] In Ref. [4], the authors investigated the properties of magnon edge states in a ferromagnetic honeycomb spin lattice with a Dzyaloshinskii–Moriya interaction (DMI). Beginning from an isotropic Heisenberg model of localized spin moments in a honeycomb lattice, Fransson et al.[5] addressed the possibility of emerging Dirac magnons. Ferreiros[7] et al. studied the influence of lattice deformations on the magnon physics of a honeycomb ferromagnet when a DMI is considered. Owerre[8] showed that the magnon Hall effect is realizable in a two-band model on the honeycomb lattice, and studied the thermal Hall conductivity. Employing the Landau–Lifshitz–Gilbert phenomenology, bulk and edge spin transport in topological magnon insulators have been studied.[9] The spin Nernst effect of magnons in a honeycomb antiferromagnet in the presence of a DMI has been demonstrated.[10,11]

Another direction that the study of topological phases has taken in recent years is manipulating the topological and transport properties of quantum systems under the influence of a periodic drive.[1517] Changes in the band structure from a nontopological band structure to a topological one can occur. In cold atom systems, periodic changes in the laser fields establish an optical lattice potential.[18,19] In solid state systems, periodically modulated quantum systems can be effectively described by a static Hamiltonian. Therefore, magnon systems deserve consideration. In this work, we study the magnon on a honeycomb lattice with a DMI under the irradiation of light.[2024] We study the effects of the interplay between the DMI and light.

The rest of this paper is organized as follows. In Section 2, we present the tight-binding magnon model on the honeycomb lattice. Based on the Brillouin–Wigner theory, we obtain the effective tight-binding Hamiltonian with a circularly polarized light. In Section 3, the Berry curvature and Chern number of the system are presented. In Section 3, the corresponding low-energy Dirac model is studied. Our conclusions are given in Section 5.

2. Model

We consider a ferromagnetic model on a honeycomb lattice. The corresponding Hamiltonian is where is the spin moment at site i. is the nearest neighbor ferromagnetic exchange coupling. The second term is the DMI between the second nearest neighbors. The constants with the vectors connecting site i to its second nearest neighbor j is shown in Fig. 1.

Fig. 1. The honeycomb lattice structure, with the nearest neighbor and the second nearest neighbors labeled by ai and di, respectively.

Using Holstein–Primakoff transformation , , and with , we obtain where ci and are magnon operators, and label the two sublattices of the honeycomb lattice denoted by different colors in Fig. 1. The Hamiltonian in the momentum space can be written as , where and where and where , , , , and . We have set the lattice constant as 1 throughout the paper. In the absence of DMI, the magnon bands are massless at the two Dirac points and .

Next, we consider the effect of the circularly polarized light represented by the time-dependent vector potential where , and E0 and ω are the field strength and frequency, respectively. Then, the response of the light-illuminated system can be obtained by substituting with . Based on the Brillouin–Wigner theory,[25,26] we can write the effective Hamiltonian as where . The second term represents the process of including the emission and absorption of the n-photon in the off-resonant light. In the off-resonant condition where the light does not cause any transition, one photon absorption and emission are enough for virtual process. Straightforward calculations show us that where is equal to the bare Hamiltonian in the absence of light in Eq. (3) except for the constant term with renormalization of the parameters: and . , and denotes the n-th Bessel function.

In addition, where , , and , and

In Fig. 2(a), we find that finite DMI opens a gap, thereby endowing the magnon bands with a nontrivial topology. Based on Eq. (8), we see that the effect of circularly polarized light can result in a gapped system. The interplay between the DMI and light can induce closing and reopening of the band gap. In Fig. 2(b), we show that the critical points where the interplay makes the energy gap are closed.

Fig. 2. Bulk magnon bands for momenta in the y direction for D=0.1 J when (a) and (b) .
3. Berry curvature and Chern number

The Berry curvature for the n-th magnon band can be calculated from The summation is over all occupied bands below the bulk gap, and is the velocity operator along the direction. En is the eigenvalue of the system (Eq. (8)). In Fig. 3, we show the Berry curvature of the lower band at D=0.1J for the varying light parameter . We can see that the Berry curvatures change with . We find that the Berry curvature increases as increases for while the reverse occurs for . The associated Chern number is defined as where the integral is over the first Brillouin zone. Consequently, the system changes from one topological magnon insulator with Chern number 1 to another one with Chern number −1. In other words, a light intensity field can lead to a topological phase transition in a topological magnon insulator.

Fig. 3. Berry curvature of the lower magnon band for D=0.1J.
4. Low-energy Dirac theory

To better understand the topological properties, we present a theoretical analysis based on an effective Hamiltonian. Expanding this Hamiltonian around the Dirac points to , we obtain the effective Dirac Hamiltonian as where corresponds to the two Dirac points, are the Pauli matrices of the sublattice pseudospin, , and . The system can be described by a static Hamiltonian in the high-frequency regime.[15,27] Then, the static time-independent Hamiltonian is written as The band dispersion relationship of the above Hamiltonian is and the the corresponding wave functions are where , , and . Based on Eq. (14), we obtain One sees that an extra light-induced term is introduced, and its sign depends on the polarization of light. This new term is clearly different from the intrinsic DMI term, and it can be defined as the light-induced effective DMI. The overall value of the DMI will decrease (increase) if the DMI term and the polarization of light have the same (different) sign. For the same sign, the critical point is defined at where the overall DMI vanishes and the band gap closes. For example, when increasing the circularly polarized light intensity, Figure 2 shows a band gap closing. The Dirac theory is consistent with the results from tight-binding theory.

The direct correspondence between the Chern number and the Hall conductance for 2D system is characterized by . The band gap ( ) is closed and reopened for fixed when varying A, and the direction of the Hall current is changed. Therefore, the direction of the Hall current can be controlled by illuminating with light due to the fact that they lead to the change in sign of the Berry curvature of the valence band. It should be noted that in our paper, we only consider . When , . The band will not close and reopen. The topological phase transition and change in direction of the Hall current will not occur.

5. Conclusion and summary

To conclude, we analyzed the topological properties of magnon in an insulator on a honeycomb lattice with DMI under the application of a circularly polarized light. We employed Brillouin–Wigner theory to obtain the tight-binding model with a circularly polarized light at the high-frequency regime. We found that the light can induce the closing and reopening of the band gap. The study about the corresponding Berry curvature and Chern number confirmed the light-induced topological phase transition. We obtained the low-energy Dirac model to better know the topological properties. The topological phase transition was light-induced handedness-dependent. The effective DMI and intrinsic DMI were responsible for the topological phase transition.

Reference
[1] Hasan M Z Kane C L 2010 Rev. Mod. Phys. 82 3045
[2] Qi X L Zhang S C 2011 Rev. Mod. Phys. 83 1057
[3] Wang X S Su Y Wang X R 2017 Phys. Rev. 95 014435
[4] Pantaleón P A Xian Y 2017 J. Phys.: Condens. Matter 29 295701
[5] Fransson J Black-Schaffer A M Balatsky A V 2016 Phys. Rev. 94 075401
[6] Kim S K Ochoa H Zarzuela R Tserkovnyak Y 2016 Phys. Rev. Lett. 117 227201
[7] Ferreiros Y Vozmediano M A H 2018 Phys. Rev. 97 054404
[8] Owerre S A 2016 J. Appl. Phys. 120 043903
[9] Rückriegel A Brataas Duine R A 2018 Phys. Rev. 97 081106
[10] Cheng R Okamoto S Xiao D 2016 Phys. Rev. Lett. 117 217202
[11] Zyuzin V A Kovalev A A 2016 Phys. Rev. Lett. 117 217203
[12] Zhang L F Ren J Wang J S Li B 2013 Phys. Rev. 87 144101
[13] Castro Neto A H Guinea F Peres N M R Novoselov K S Geim A K 2009 Rev. Mod. Phys. 81 109
[14] Kane C L Mele E J 2005 Phys. Rev. Lett. 95 146802
[15] Lindner N H Rafael G Galitski V 2011 Nat. Phys. 7 490
[16] Kitagawa T Berg E Rudner M Demler E 2010 Phys. Rev. 82 235114
[17] Rudner M S Lindner N H Berg E Levin M 2013 Phys. Rev. 3 031005
[18] Jotzu G Messer M Desbuquois R Lebrat M Uehlinger T Greif D Esslinger T 2014 Nature 515 237
[19] Bilitewski T Cooper N R 2015 Phys. Rev. 91 063611
[20] Zhong W T Liu J Hua D Hou H Y Yan K J 2018 Acta Phys. Sin. 18 184208 (in Chinese)
[21] Liu M J Zhang X W Wang J Qin Y B Chen Y Hua Huang W 2018 Acta Phys. Sin. 20 207801 (in Chinese)
[22] Wang X Pu H B Liu Q An L Q 2018 Chin. Phys. 27 108502
[23] Qi K Z Tang Y Yang L P Yang F Y Guo W J 2018 Chin. Phys. 27 107801
[24] Shen N C 2016 Physics 45 790 (in Chinese)
[25] Mikami T Kitamura S Yasuda K Tsuji N Oka T Aoki H 2016 Phys. Rev. 93 144307
[26] Mohan P Saxena R Kundu A Rao S 2016 Phys. Rev. 94 235419
[27] Oka T Aoki H 2009 Phys. Rev. 79 081406